Enstrophy and dissipation must have the same scahng exponent in the high Reynolds 

number limit of fluid turbulence 
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Writing tlie Poisson equation for the pressure in the vorticity-strain form, we show that the pres- 
sure has a finite inertial range spectrum in the high Reynolds number limit of isotropic turbulence 
only if the anomalous scaling exponents fi and ^uj for the dissipation and enstrophy (squared vortic- 
ity) are equal. Since a finite inertial range pressure spectrum requires only very weak assumptions 
about high Reynolds number turbulence, we conclude that the inference from experiment and direct 
numerical simulation that these exponents are different must be a finite range scaling result which 
will not survive taking the high Reynolds number limit. 
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It has long been recognized ^ that the local rate of 
dissipation in high Reynolds number turbulence exhibits 
anomalous scaling, characterized by the fluctuations in 
er{x,t), the locally averaged dissipation rate over a re- 
gion of size r, where e is the local rate of dissipation per 
unit mass of turbulent kinetic energy. Here we find it 
more convenient to work with the spectrum of dissipation 
fluctuations in the inertial range, which is conventionally 
assumed to have the form 



(1) 



where fi is the putatively universal exponent describing 
the dissipation fluctuations in high Reynolds number tur- 
bulence [||,(e) is the average rate of energy dissipation 
per unit mass, k is the wave number, L is the geomet- 
rically defined external length scale, and Ce is a dimen- 
sionless constant which is probably not universal. 

For later purposes, note that the dissipation rate e is 
related to the symmetric part of the strain tensor as 



where = SijSji, Sij — l/2(dui/ dxj + duj/dxi), and 
V is the kinematic viscosity of the fluid. The relationship 
between S*, the antisymmetric part of the strain tensor 
Gij = l/2{dui/dxj —duj/dxi) and the enstrophy fl = u'^, 



n 

2 



dui duj 

dxj dxi 



(2) 



Using the Navier-Stokes equations, (g) can also be related 
to the pressure p via 



or 



^2 duj 
^ dxj dxi ' 



(3) 



(4) 



and the density of the fluid is taken as p = 1 with no loss 
of generality. The derivation of (||), (||) and ^ makes 
essential use of incompressibility. 

For the enstrophy, a quantity analagous to the dissi- 
pation can be defined as 



D = vuj\ 



(5) 



The anomalous scaling of D is exhibited in its power spec- 
trum, which is assumed to have the inertial range form 



ED{k)^CD{e?k-\kLy 



(6) 



where pi^ is a scaling exponent describing the anomalous 
scaling of enstrophy fluctuations. 

There has been an interesting recent controversy con- 
cerning the relative scaling of enstrophy and dissipation. 
Older direct numerical simulations (DNS) ||^ showed 
that vorticity was more intermittent than strain for single 
point quantities. The one-dimensional measurements of 
the streamwise components of e and w, obtained at both 
high and low Reynolds numbers ||^, conclude that the 
degree of intermittency in the dissipation and enstrophy 
fields are not the same. On the other hand L'vov and 
Procaccia have argued on general symmetry grounds 
that the asymptotic scaling exponents must be equal; 
i.e. that = /i. This motivated Chen, Sreenivasan and 
Nelkin |^ to study the inertial range scaling of enstrophy 
and dissipation in a DNS at moderate Reynolds number 
R\ ~ 216. They found that the enstrophy was clearly 
more intermittent than the dissipation suggesting that 
/icj > /i. Finally, He et al |^ have studied the statistics 
of dissipation and enstrophy induced by a set of Burgers 
vortices. They found that finite-range scaling exponents 
are different for these two quantities, but that for this 
model system, the asymptotic scaling exponents for dis- 
sipation and enstrophy can be shown to be equal in the 
limit of infinite Reynolds number. It is then of consider- 
able interest to see if this can be shown more generally 
without resort to a particular model system. 
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In this paper I show on quite general grounds, starting 
from that /Xoj = /Lt. I use only the assumption that 
the pressure has a finite power spectrum in the inertial 
range when v ^ 0. The qualitative idea is that each of 
the source terms on the right hand side of (^) scales as 
(l/i^) but their difference must be independent of i^. This 
can only be true if each source term separately has the 
same scaling exponent. 

In writing (|l|) and the idea of an inertial range 
has been mentioned, but it is important to define it more 
precisely. The basic assumption is that certain quan- 
tites are independent of the kinematic viscosity v in the 
limit that — > 0. This is conventionally assumed for 
the average dissipation rate (e) , and the usual picture of 
dissipation fluctuations assumes that it is also true for 
the spectrum of dissipation flucutations. This is indi- 
cated in (0) and (|^) by the appearance of the extenal 
length scale L rather than the dissipation length scale 
77. It is conventionally assumed that the pressure spec- 
trum has a finite inertial range. This follows naturally 
from the Poisson equation in the form Hill and 

Wilczak have shown that the pressure structure func- 
tion {[p{x + r) — can be written as an integral 
over the three independent fourth order velocity struc- 
ture functions {[u{x + r) ~ u{x)]'^), {[v{x + r) ~ v{x)]'^), 
and {[u{x + r) — u{x)]'^[v{x + r) — v{x)]'^), where u{x) 
is the longitudinal component of the velocity along the 
line joining the points x and (x + r), and v{x) is one of 
the transverse components of velocity perpindicular to 
this line. If these three structure functions have finite 
inertial ranges scaling approximately in the Kolmogorov 
form r^/^, then the integrals in Hill and Wilczak are com- 
fortably convergent, and the pressure structure function 
and pressure spectrum have finite inertial ranges. Hill 



Ll[ have recently 
ate the pressure 



and Boratav and Nelkin and Chen 
used the Hill- Wilczak formula to calcu 
structure function for DNS and for experiments. The 
resulting expression exhibits interesting partial cancel- 
lations, and a high degree of sensitivity to the detailed 
form of the fourth order velocity structure functions, but 
there are no singular contributions and no mathematical 
difficulties. 

Starting from (^), however, great care must be taken 
to get a finite pressure spectrum in the limit 1/ — > 0. 
Taking the fourier transform of (^, taking the absolute 
square, and using (|^) and (^), the pressure spectrum can 
be written as 



:[CAkL)-^ + CoikL)-^- - F,{kL)], (7) 



where the last term on the right hand side of (0) repre- 
sents the cross term between enstrophy and dissipation. 
The most natural assumption is that the cross term scales 
as 



F4kL) = C^kLY 



(8) 



scaling dimension and fi/2 respectively, in which 

case the exponent fi^ is given by 



(9) 



If the pressure spectrum is to be finite in the limit 
v the right hand side of (|7|), which is proportional 

to v^^ must cancel identically for all values of k. If the 
cross term satisfies (^, this is only possible if the scaling 
exponents in (|^) satisfy 



(10) 



From a field theoretical viewpoint, it is natural to think 
of the enstrophy and dissipation as scaling variables with 



so that the scaling exponents for enstrophy and dissipa- 
tion must be equal. 

If the cross term has a more complicated structure, and 
does not satisfy a power law, then the argument given 
above no longer has a simple algebraic structure. Taking 
the Fourier transform of (ra), however, it is very hard to 
see how a finite value of the inertial range pressure spec- 
trum can arise without an almost total cancellation of 
the two source terms on the right hand side. Such a can- 
cellation seems implausible unless the two source terms 
have the same scaling exponent. 

At this stage of the argument, this cancellation of 
the dominant terms gives a pressure spectrum which is 
identically zero. This just emphasizes that (||) is not 
a good starting point for calculating the pressure spec- 
trum. In fact, both (|^) and (||) also have non-singular 
background terms which have been omitted here. In the 
quasi-Gaussian approximation, these background terms 
will scale as u^t^/^k^/^ ^ and will give a non- vanishing 
contribution to the pressure spectrum of order e^/^fc~^/^ 
as expected from conventional dimensional arguments of 
the 1941 Kolmogorov type In the individual spectra 
(^, and (|^), however, these background terms are neg- 
ligible in the inertial range compared to the dominant 
anomalous scaling terms. It is only after the dominant 
singularity has been completely cancelled in the pressure 
that these remaining terms become important. As dis- 
cussed in pO| and pT| ] these remaining terms are sensitive 
to deviations from the quasi-Gaussian approximation so 
that the pressure spectrum and structure function are 
able to probe detailed properties of the fourth order ve- 
locity structure frunctions. 

Finally I would add two statements about the underly- 
ing physics. First I would emphasize that the result pre- 
sented here applies only in the limit of very high Reynolds 
number. As discussed in the apparent scaling can 
be different for finite but very large Reynolds number. 
When discussing subtle questions of small differences in 
scaling, experiment may be far from asymptotic, even 
for atmospheric turbulence. Second, I find it interesting 
that the basic result derived here from simple scaling ar- 
guments reminiscent of statistical field theory has been 
also been derived in Q from considerations of localized 
vortex structures. It is rare in turbulence theory that the 
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same result is obtained from the statistical point of view 
and from the point of view of vortex structures. 

The idea for this short derivation occurred at the meet- 
ing "Turbulence: Challenges for the 21st Century," at 
Los Alamos, May 18-21, 1998. This meeting was also in 
celebration of Robert H. Kraichnan's 70th birthday. A 
key suggestion was made by Boris Shraiman, who pointed 
out to me at that meeting that the right hand side of 
must cancel almost everywhere. From this I have 
deduced, with quite modest assumptions, that the dif- 
ferences in scaling of enstrophy and dissipation observed 
by various authors must be a finite range scaling effect. 
This is in agreement with the recent model calculation of 
He et al ^ , and can be thought of as a modest birthday 
gift from myself to Bob Kraichnan. 
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